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^O ■ QCD-based thermodynamics at zero and finite quark chemical potential is 

studied using an extended Nambu and Jona-Lasinio approach in which quarks 
couple simultaneously to the chiral condensate and to a background temporal 
gauge field representing Polyakov loop dynamics. This so-called PNJL model 
thus includes features of both deconfinement and chiral symmetry restoration. 
We discuss the phase diagram as it emerges from this approach in close com- 
parison with results from lattice QCD thermodynamics. The critical point, 
separating crossover from first order phase transition, is investigated with 
special focus on its quark mass dependence, starting from the relatively large 
masses presently accessible by lattice simulations, down to the chiral limit. 

1 Introduction and Basics 

Models of the Nambu and Jona-Lasinio (N JL) type [1] have a long history and have 
been used extensively to describe the dynamics and thermodynamics of the lightest 
hadrons [2-5], including investigations of phase diagrams [6,7]. Such schematic 
models offer a simple and practical illustration of the basic mechanisms that drive 
spontaneous chiral symmetry breaking, a key feature of QCD in its low-temperature, 
low-density phase. 

The NJL model is based on an effective Lagrangian of relativistic fermions 
(quarks) which interact through local current- current couplings, assuming that glu- 
onic degrees of freedom can be frozen into pointlike effective interactions between 
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quarks. Lattice QCD results for the gluonic field strength correlation function [9] 
demonstrate that the colour correlation length, i.e. the distance over which colour 
fields propagate in the QCD vacuum, is small, of order 0.2 fm corresponding to a 
characteristic momentum scale A of order 1 GeV. Consider now the basic non-local 
interaction between two quark colour currents, J? = ip'-fHiip , where ti are the gen- 
erators of the SU(N C ) colour gauge group. The contribution of this current-current 
coupling to the action is: 

S mt = ~Jd A xd 4 yJt(x)g 2 D%(x,y)J»(y) , (1) 

where D l j v is the full gluon propagator and g is the QCD coupling. In perturbative 
QCD this Si n t generates the familiar one-gluon exchange interaction between quarks 
and maintains its non-local structure. That is the situation realised in the quark- 
gluon phase at extremely high temperatures. As one approaches the hadronic phase 
around a critical temperature of about 0.2 GeV, the propagating gluons experience 
strong screening effects which cannot be handled perturbatively. If the range over 
which colour can be transported is now restricted to the short distance scale A^ 1 , 
while typical momentum scales (Fermi momenta) of the quarks are small compared 
to A, then the quarks experience an interaction which can be approximated by a 
local coupling between their colour currents: 

C mt = -G c j;(x)J?(x) , (2) 

where G c ~ g 2 Ar 2 is an effective coupling strength of dimension length 2 which en- 
codes the QCD coupling, averaged over the relevant distance scales, in combination 
with the squared correlation length, A -2 . In essence, by "integrating out" gluon de- 
grees of freedom and absorbing them in the four-fermion interaction Ci nt , the local 
SU(N C ) gauge symmetry of QCD is now replaced by a global SU(N C ) symmetry 
of the NJL model. Apart from this step, the interaction Lagrangian (2) evidently 
preserves the chiral SU(Nf) x SU(Nf) symmetry that it shares with the original 
QCD Lagrangian for Nf massless quark flavours. 

A Fierz transform of the colour current-current interaction (2) produces a set of 
exchange terms acting in quark-antiquark channels. For the Nf = 2 case: 

£n*~>f [(W 2 + (#76*0) 2 ]+-, ( 3 ) 

where f = (ti,T2,Ts) are the isospin SU(2) Pauli matrices. Not shown for brevity 
is a series of terms with combinations of vector and axial vector currents, both in 
colour singlet and colour octet channels. The constant G is proportional to the 
colour coupling strength G c . Their ratio is uniquely determined by N c and Nf. 



Eq.(3) is the starting point of the standard NJL model. In the mean field 
(Hartree) approximation, the NJL equation of motion leads to the gap equation 

m = m — G(ijjip) . (4) 

With a small bare (current) quark mass itlq as input, this equation links the dynam- 
ical generation of a large constituent quark mass m to spontaneous chiral symmetry 
breaking and the appearance of the quark condensate 

<*» = -1> Jjm <7V<0)**)> = -™,* J gj TjZ^ ■ (5) 

For mo = a non-zero quasiparticle mass develops dynamically, together with a 
non-vanishing chiral condensate, once G exceeds a critical value. The procedure 
requires a momentum cutoff A ~ 2m beyond which the interaction is "turned off" . 
Note that the strong interaction, by polarizing the vacuum and turning it into a 
condensate of quark-antiquark pairs, transforms an initially pointlike quark with its 
small bare mass itlq into a massive quasiparticle with a finite size. (Such an NJL- 
type mechanism is commonly thought to be at the origin of the phenomenological 
constituent quark masses m ~ 0.3-0.4 GeV). 

While the NJL model illustrates the transmutation of originally light (or even 
massless) quarks and antiquarks into massive quasiparticles, it generates at the same 
time pions as Goldstone bosons of spontaneously broken chiral symmetry. NJL type 
approaches have also been used extensively to explore colour superconducting phases 
at high densities through the formation of various sorts of diquark condensates [6, 8]. 

Despite their widespread use, NJL models have a principal deficiency. The reduc- 
tion to global (rather than local) colour symmetry has the consequence that quark 
confinement is missing. Confinement is the second key feature of low-energy QCD 
besides spontaneous chiral symmetry breaking. While confinement is a less signif- 
icant aspect for N c = 2 thermodynamics which can be described quite successfully 
using the simplest NJL approach [10], it figures prominently for N c = 3 QCD. There 
have been ongoing discussions whether deconfinement and the restoration of chiral 
symmetry are directly connected in the sense that they appear at the same transi- 
tion temperature T c , as suggested by earlier lattice computations. In any case, as 
one approaches T c from above, all versions of the "classic" NJL model encounter the 
problem that they operate with the "wrong" degrees of freedom. Quarks as coloured 
quasiparticles are incorrectly permitted to propagate over large distances even in the 
hadronic sector of the phase diagram. In contrast, confinement and spontaneous chi- 
ral symmetry breaking imply that QCD below T c turns into a low-energy effective 
theory of weakly interacting Goldstone bosons (pions) with derivative couplings to 
colour-singlet hadrons (rather than quarks). 

In the limit of infinitely heavy quarks, the deconfinement phase transition is 
characterized by spontaneous breaking of the Z(3) center symmetry of QCD. The 



corresponding order parameter is the thermal Wilson line, or Polyakov loop, winding 
around the imaginary time direction with periodic boundary conditions: 
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with j3 = 1/T the inverse temperature. Here A 4 = iA° is the temporal component 
of the Euclidean gauge field (A, A 4 ) and V denotes path ordering. In the presence 
of dynamical quarks the Z(3) symmetry is explicitly broken. The Polyakov loop 
ceases to be a rigourous order parameter but still serves as an indicator of a rapid 
crossover towards deconfinement. 

Recent developments have aimed at a synthesis of the NJL model with Polyakov 
loop dynamics. The principal idea is to introduce both the chiral condensate (ipip) 
and the Polaykov loop (f> as classical, homogeneous fields which couple to the quarks 
according to rules dictated by the symmetries and symmetry breaking patterns of 
QCD, thus unifying aspects of confinement and chiral symmetry breaking. We 
refer to this combined scheme as the PNJL (Polyakov-loop-extended NJL) model. 
The present writeup is largely based on our previous ref.[ll]. It reviews results of 
calculations as well as presenting an outlook for further steps to be pursued in close 
comparison with results from lattice QCD thermodynamics. 

2 Introducing the PNJL model 

Throughout this presentation we work with two flavours (Nf = 2) and specify the 
PNJL Lagrangian [11] as follows. Its basic ingredients are the Nambu and Jona- 
Lasinio type four-fermion contact term and the coupling to a (spatially constant) 
temporal background gauge field representing Polyakov loop dynamics: 

Unjl = 4> (i-Y^D" - m ) ^ + ^ [(^) 2 + (^7 5 rV>) 2 ] - U (</>[A], cf>*[A}; T) , (7) 

where if) = (if) u , ipa) is the quark field, 

D^ = d^-iA^ and A* = 8$A° , (8) 

with A = —iA 4 . The gauge coupling g is conveniently absorbed in the definition 
of A^(x) = gA^(x)^f where A% is the SU(3) gauge field and A a are the Gell-Mann 
matrices. The two-flavour current quark mass matrix is itiq = diag(m u ,md) and 
we shall work in the isospin symmetric limit with m u = rrid = itlq. As previously 
mentioned, G is the coupling strength of the chirally symmetric four-fermion inter- 
action. 



The effective potential W(0, 0*;T) is expressed in terms of the traced Polyakov 
loop (6), reduced to our case of a constant Euclidean field A 4 : 



-Tr c exp 



iA 4 



(9) 



In a convenient gauge (the so-called Polyakov gauge), the Polyakov loop matrix can 
be given a diagonal representation [12]. 

The effective potential U has the following general features. It must satisfy the 
Z{3) center symmetry just like the pure gauge QCD Lagrangian. Furthermore, in 
accordance with lattice results for the behaviour of the Polyakov loop as a function 
of temperature T, the potential U must have a single minimum at = at small T, 
while at high T it develops a second minimum which becomes the absolute minimum 
above a critical temperature To. In the limit T — ► oo we have — »■ 1. The following 
general form is chosen for U, including a 3 term which reflects the underlying Z(3) 
symmetry: 

W(0,0*;T) _ b 2 (T) ,„, &3,, 3 , ,„- , b 4 



+ <T) + i (W (io) 



with 



T 4 2 r r 6 vr r 7 4 



6 2 (T) = a + a 1 (^]+a 2 ^] + a 3 (f ) • (H) 



A precision fit of the coefficients a^, 6j is performed to reproduce the pure-gauge 
lattice data. 

There is a subtlety about the Polyakov loop field, 0, and its conjugate, 0*, in the 
presence of quarks. At zero chemical potential we have = 0*, i.e. the field is real, 
it serves as an order parameter for deconfinement and a mean-field calculation is 
straightforward. At non-zero quark chemical potential, Z(3) symmetry is explicitly 
broken and differs from 0* while their thermal expectation values (0) and (0*) 
remain real [13] . A detailed analysis of the stationary points of the action under these 
conditions requires calculations beyond mean field which will be reported elsewhere 
[14]. We proceed here, as in [11], by introducing $ = (0) and $ = (0*) as new 
independent field variables which replace and 0* in Eq.(10). This approximate 
prescription corresponds to a modified mean-field scheme which can account for the 
difference between $ and $ in the presence of quarks. The more accurate treatment 
is under way. 

Using standard bosonization techniques, we introduce the auxiliary bosonic fields 
a and -jf for the scalar-isoscalar and pseudoscalar-isovector quark bilinears in Eq.(7). 
The expectation value of the a field is directly related to the chiral condensate by 
(a) = G(ipip) and the gap equation becomes 

m = m — (a) . (12) 



a 


a\ 


a 2 


a 3 


h 


h 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



Table 1: Parameter set used in [11] for the Polyakov loop potential (10, 11). 



Note that (a) is negative in our representation, and the chiral (quark) condensate 

is (-0V) = (ii + VW^)- 

Before passing to the actual calculations, we summarize basic assumptions be- 
hind Eq.(7) and comment on limitations to be kept in mind. The PNJL model 
reduces gluon dynamics to a) chiral point couplings between quarks, and b) a sim- 
ple static background field representing the Polyakov loop. This picture can be 
expected to work only within a limited range of temperatures. At large T, trans- 
verse gluons are known to be thermodynamically active degrees of freedom, but they 
are ignored in the PNJL model. To what extent this model can reproduce lattice 
QCD thermodynamics is nonetheless a relevant question. We can assume that its 
range of applicability is, roughly, T < (2 — 3)T C , based on the conclusion drawn in 
ref. [15] that transverse gluons start to contribute significantly for T > 2.5 T c . 

3 Parameter fixing 

The parameters of the Polyakov loop potential U are fitted to reproduce the lat- 
tice data [16] for QCD thermodynamics in the pure gauge sector. Minimizing 
W($,$,T) one has $ = $ and the pressure of the pure-gauge system is evalu- 
ated as p(T) = —U(T) with $(T) determined at the minimum. The entropy and 
energy density are then obtained by means of the standard thermodynamic relations. 
Fig. 1(a) shows the behaviour of the Polyakov loop as a function of temperature, while 
Fig. 1(b) displays the corresponding (scaled) pressure, energy density and entropy 
density. The lattice data are reproduced extremely well using the ansatz (10,11), 
with parameters summarized in Table 1. The critical temperature To for decon- 
finement appearing in Eq. (11) is fixed at T = 270 MeV in the pure gauge sector. 



The pure NJL model part of the Lagrangian (7) has the following parameters: 
the "bare" quark mass m , a three-momentum cutoff A and the coupling strength 
G. We fix them by reproducing the known chiral physics in the hadronic sector at 
T = 0: the pion decay constant /„-, the chiral condensate K^puipu)] 1 ^ 3 and the pion 
mass m n are evaluated in the model and adjusted at their empirical values. The 
results are shown in Table 2. 
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Figure 1: (a): Polyakov loop as a function of temperature in the pure gauge sector, compared to 
corresponding lattice results taken from Rcf. [17]. (b): Scaled pressure, entropy density and energy 
density as functions of the temperature in the pure gauge sector, compared to the corresponding 
lattice data taken from Rcf. [16]. 



A [GeV] 


G[GeV- 2 ] 


m [MeV] 


0.651 


10.08 


5.5 


|(VvM| 1/3 [MeV] 


U[MeV] 


77V [Me V] 


251 


92.3 


139.3 



Table 2: Parameter set used for the NJL model part of the effective Lagrangian (7), 
and the resulting physical quantities. These values of the parameters yield a con- 
stituent quark mass m = 325 MeV. 



4 Finite \i thermodynamics 

4.1 General features 

We now extend the model to finite temperature and chemical potentials using the 
Matsubara formalism. We consider the isospin symmetric case, with an equal num- 
ber of u and d quarks (and therefore a single quark chemical potential //). The 
quantity to be minimized at finite temperature is the thermodynamic potential per 
unit volume: 

2 

n = w($,$,t) + ^ 

_ 2N f T I ^ {Tr c In [l + e -(^-£)/ r ] + Tr c In [l + e^ E "^/ T ] } 

- W f Jj^E p 6{tf-f) , (13) 



where \x — // + iA 4 and E p = yp 2 + m 2 is the quark quasiparticle energy. The last 
term involves the NJL three- momentum cutoff A. The second (finite) term does not 
require any cutoff. 

Notice that the coupling of the Polyakov loop to quarks effectively reduces 
the residues at the quark quasiparticle poles as the critical temperature is ap- 
proached from T > T c : expanding the logarithms in the second line of (13) one 
finds Tr c ln (1 + exp[— (E p — \x — iA±)/T\) = 30exp[— (E p — fx)/T\ + ... , with <fi 
then to be replaced by (0) = $ which tends to zero as T — > T c . 

From the thermodynamic potential (13) the equations of motion for the mean 
fields a, $ and $ are determined through 

dtt dtt dtt 

a? = - a¥ = ' »=°- < 14 » 

This set of coupled equations is then solved for the fields as functions of temperature 
T and quark chemical potential /x. Fig. 2(a) shows the chiral condensate together 
with the Polyakov loop $ as functions of temperature at \x — where we find $ = $. 
One observes that the introduction of quarks coupled to the a and $ fields turns 
the first-order transition seen in pure-gauge lattice QCD into a continuous crossover. 
The crossover transitions for the chiral condensate (ipip) and for the Polyakov loop 
almost coincide at a critical temperature T c ~ 220 MeV (see Fig. 2(b)). We point 
out that this feature is obtained without changing a single parameter with respect 
to the pure gauge case. The value of the critical temperature found here is a little 
high if compared to the available data for two-flavour Lattice QCD [18] which give 
T c = (173 ± 8) MeV. For quantitative comparison with existing lattice results we 
choose to reduce T c by rescaling the parameter To from 270 to 190 MeV. In this 
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Figure 2: Left: scaled chiral condensate and Polyakov loop $(T) as functions of 
temperature at zero chemical potential. Right: plots of d(ipif})/dT and dQ/dT. 

case we loose the perfect coincidence of the chiral and deconfinement transitions, 
but they are shifted relative to each other by less than 20 MeV. When defining T c in 
this case as the average of the two transition temperatures we find T c = 180 MeV. 

4.2 Detailed comparison with lattice data 

The primary aim is now to compare predictions of our PNJL model with the lattice 
data available for full QCD thermodynamics (with quarks included) at zero and 
finite chemical potential ji. Consider first the pressure p(T, fj, — 0) = —Q (T, \x = 0) 
of the quark-gluon system at zero chemical potential. Our results are presented in 
Fig. 3 in comparison with corresponding lattice data. We point out that the input 
parameters of the PNJL model have been fixed independently in the pure gauge 
and hadronic sectors, so that the calculated pressure is a prediction of the model, 
without any further tuning of parameters. With this in mind, the agreement with 
lattice results is quite satisfactory. Also shown in Fig. 3 is the result obtained in the 
standard NJL model. Its deficiencies are evident. At low temperatures the pressure 
comes out incorrect. The missing confinement permits quarks to be active degrees 
of freedom even in the forbidden region T < T c . At high temperatures, the standard 
NJL result for the pressure is significantly lower than the one seen in the lattice 
data. The gluonic thermodynamics is missing altogether in the NJL model, whereas 
in the PNJL model it is partially taken into account by means of the Polyakov loop 
effective potential W($,$,T). As stated previously, the range of validity of this 
approach is limited, however, to temperatures smaller than 2.5 T c , beyond which 
transverse gluon degrees of freedom become important. 

The introduction of the Polyakov loop within the PNJL quasiparticle model 
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Figure 3: Scaled pressure divided by the Stefan-Boltzmann (ideal gas) limit as a function of 
temperature at zero chemical potential: comparison between our PNJL model prediction [11] (full 
line), the standard NJL model result (dashed) and lattice results corresponding to Nt = 4 and 
Nt = 6. Lattice data are taken from Ref. [19] 

leads to a remarkable improvement in basically all thermodynamic quantities. The 
coupling of the quark quasiparticles to the field $ reduces their weight as thermody- 
namically active degrees of freedom when the critical temperature T c is approached 
from above. The quasiparticle exponentials exp[— (E p ±li)/T] are progressively sup- 
pressed in the thermodynamic potential as T — > T c . This is what can be interpreted 
as the effect of confinement in the context of the PNJL model. 

One must note that the lattice data are grouped in different sets obtained on 
lattices with temporal extent N t = 4 and N t = 6, both of which are not continuum 
extrapolated. In contrast, our calculation should, strictly speaking, be compared 
to the continuum limit. In order to perform meaningful comparisons, the pressure 
is divided by its asymptotic high-temperature (Stefan-Boltzmann) limit for each 
given case. At high temperatures our predicted curve should be located closer to 
the N t — 6 set than to the one with N t = 4. This is indeed the case. 

At non-zero chemical potential, quantities of interest that have become accessible 
in lattice QCD are the "pressure difference" and the quark number density. The 
(scaled) pressure difference is defined as: 

Ap(T,ti) P {T,fjL)-p{T,fi = 0) 



rp^ 



rp4 



(15) 



A comparison of Ap, calculated in the PNJL model, with two-flavour lattice results 
is presented in Fig. 4. This figure shows the scaled pressure difference as a function 
of the temperature for a series of chemical potentials, with values ranging between 
jx = 0.2 Tc and ii ~ T c where T[ \ = T c (p = 0). The agreement between our 
results [11] and the lattice data is quite satisfactory 
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Figure 4: Scaled pressure difference as a function of temperature at different values of the quark 
chemical potential (results from Ref. [If]), compared to lattice data taken from Ref. [20]. 
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Figure 5: Scaled quark number densities [11] as a function of temperature at different values of 
the chemical potential, compared to lattice data taken from Ref. [20]. 
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A related quantity for which lattice results at finite fi exist, is the scaled quark 
number density, defined as: 

n q (T,fi) 1 dQ{T,fi) 

—T^ = -^^d^- (16) 

Results [11] for n q as a function of the temperature, for different values of the 
quark chemical potential, are shown in Fig. 5 in comparison with corresponding 
lattice data [20] . Also in this case, the agreement between our PN JL model and the 
corresponding lattice data is surprisingly good. 

It is a remarkable feature that the quark densities and the pressure difference at 
finite fi are so well reproduced even though the lattice "data" have been obtained by 
a Taylor expansion up to fourth order in fi, whereas our thermodynamical potential 
is used with its full functional dependence on fi. We have examined the convergence 
in powers of \x by expanding Eq. (13). It turns out that the Taylor expansion to 
order fi 2 deviates from the full result by less than 10 % even at a chemical potential 
as large as fi ~ T c . When expanded to 0(fi 4 ), no visible difference is left between 
the approximate and full calculations for all cases shown in Figs. 4 and 5. 

An exact copy of our PNJL model [11] has recently been employed in [21] to 
calculate susceptibilties and higher order derivatives in the expansion of the pressure 
p(T,fi) = —Q(T,fi) around fi = 0. 



p(T, fi) 



n even 

The resulting quark number susceptibility 

c 2 (T) * < '^ 



£ c «( T ) (fT • ( 17 ) 



2T 2 \dfi 



u=0 



compares well with lattice QCD computations. The higher-order coefficients c^q 
reproduce the corresponding lattice data around T c very well, but C4 as obtained 
in the PNJL calculation tends to be too large at higher temperatures. For a more 
quantitative understanding, further steps are yet necessary towards a consistent 
treatment beyond the mean-field level [14]. 

5 Phase diagram 

Lattice data for the QCD phase diagram exist up to relatively high temperatures, but 
extrapolations to non-zero chemical potential are still subject to large uncertainties. 
It is nonetheless instructive to explore the phase diagram as calculated in the PNJL 
model [11] in comparison with present lattice QCD results. In particular, questions 
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Figure 6: Phase diagram of the PNJL model for different values of the bare quark mass mo. 
Dashed lines correspond to a first order phase transition, full lines to a crossover. The band 
represents the extrapolation from lattice QCD results [22]. Also shown for orientation are the 
chemical freezeout "data" obtained through a thermal fit [23] . 

about the sensitivity of this phase diagram with respect to changes of the input 
quark masses will be addressed. This is an important issue, given the fact that most 
lattice QCD computations so far encounter technical limitations which restrict the 
input bare quark masses to relatively large values. The PNJL approach permits 
to vary the bare quark mass in a controlled way compatible with explicit chiral 
symmetry breaking in QCD. One can therefore interpolate between large quark 
masses presently accessible in lattice simulations, the physically relevant range of 
light quark masses around 5 MeV and further down to the chiral limit. 

Fig. 6 presents our two-flavour PNJL results of the phase boundaries in the (T, /x) 
plane. These calculations should still be considered as an exploratory study since 
they do not yet include explicit diquark degrees of freedom, an important ingredient 
when turning to large chemical potentials, and the Polyakov loop fields are still 
treated in an approximate mean- field framework. Some interesting tendencies are, 
however, already apparent at the present stage. 

Curves are shown for three different values of the bare quark mass. For itiq = 50 
MeV our result falls within the broad band of lattice extrapolations using an expan- 
sion in powers of the quark chemical potential. Reducing the bare quark mass toward 
physically realistic values leaves the phase diagram at small chemical potentials basi- 
cally unchanged. However, the phase boundary is shifted quite significantly to lower 
temperatures at increasing chemical potentials when ttlq is lowered. Also shown in 
the figure is the position of the critical point separating crossover from first order 
phase transition. In fact, examining the chiral condensate and the Polyakov loop as 
functions of temperature for a broad range of chemical potentials (see Fig. 7), one 
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Figure 7: Constituent quark mass (a) and Polyakov loop (6) as functions of temperature for 
different values of the chemical potential. In both figures, mo = 5.5 MeV. 

observes that there is a critical chemical potential above which these two quantities 
indicate a discontinuous jump from the confined (chirally broken) to the deconfined 
(chirally restored) phase. While a more precise location of this critical point is sub- 
ject to further refined calculations [14], the qualitative features outlined here are 
expected to remain, such as the observation that the position of the critical point 
depends sensitively on the input quark mass. 

6 Conclusions 

The PNJL approach represents a minimal synthesis of the two basic principles that 
govern QCD at low temperatures: spontaneous chiral symmetry breaking and con- 
finement. The respective order parameters (the chiral quark condensate and the 
Polyakov loop) are given the meaning of collective degrees of freedom. Quarks cou- 
ple to these collective fields according to the symmetry rules dictated by QCD itself. 
A limited set of input parameters is adjusted to reproduce lattice QCD results 
in the pure gauge sector and pion properties in the hadron sector. Then the quark- 
gluon thermodynamics above T c up to about twice the critical temperature is well 
reproduced, including quark densities up to chemical potentials of about 0.2 GeV. 
In particular, the PNJL model correctly describes the step from the first-order de- 
confinement transition observed in pure-gauge lattice QCD (with T c ~ 270 MeV) 
to the crossover transition (with T c around 200 MeV) when Nf = 2 light quark 
flavours are added. The non-trivial result is that the crossovers for chiral symmetry 
restoration and deconfinement almost coincide, as found in lattice simulations. The 
model also reproduces the quark number densities and pressure difference at vari- 
ous chemical potentials surprisingly well when confronted with corresponding lattice 
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data. Considering that the lattice results have been found by a Taylor expansion in 
powers of the chemical potential, this agreement indicates rapid convergence of the 
power series in \x. 

The phase diagram predicted in this model has interesting implications. Starting 
from large quark masses an extrapolation to realistic small quark masses can be 
performed. The location of the critical point turns out to be sensitive to the input 
value of the bare (current) quark mass. 

The conclusion to be drawn at this point is as follows. A quasiparticle approach, 
with its dynamics rooted in spontaneous chiral symmetry breaking and confine- 
ment and with parameters controlled by a few known properties of the gluonic and 
hadronic sectors of the QCD phase diagram, can account for essential observations 
from two-flavour N c = 3 lattice QCD thermodynamics up to about twice the critical 
temperature of about 0.2 GeV. Presently ongoing further developments include: 

• systematic steps beyond the mean-field approximation; 

• extensions to 2+1 flavours; 

• inclusion of explicit diquark degrees of freedom and investigations of colour 
superconductivity in the high density domain; 

• detailed evaluations of susceptibilities and transport properties at finite chem- 
ical potential. 
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